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<^ \ We perform a comparison of the radiative and collisional parton energy losses in an ex- 
\q ; panding quark-gluon plasma. The radiative energy loss is calculated within the light-cone 
path integral approach [4]. The collisional energy loss is calculated using the Bjorken 
method with an accurate treatment of the binary collision kinematics. Our numerical 
results demonstrate that for RHIC and LHC conditions the collisional energy loss is rel- 
Ph| atively small in comparison to the radiative one. We find an enhancement of the heavy 
quark radiative energy loss as compared to that of the light quarks at high energies. 

1. The suppression of high-j>r hadrons in v4v4-collisions (usually called jet quenching) 
y—{ \ observed in the experiments at RHIC (for a review, see [1]) is widely believed to be due 
to the parton energy loss in the hot quark-gluon plasma (QGP) produced at the initial 
stage of nucleus-nucleus collision. The parton energy loss may come from the collisional 
energy loss and the induced gluon radiation. The first estimate of the collisional energy 
loss in the QGP has been done by Bjorken [2]. The radiative energy loss has been under 
active investigation in the last years [3, 4, 5, 6, 7, 8, 9] (for a review, see [10]). The 
calculations of the radiative energy loss within the BDMPS [3, 10] and the light-cone 
^ ■ path integral (LCPI) [4, 6, 7]) approaches demonstrate that for high energy partons the 
energy loss is likely dominated by the induced gluon radiation. The estimates given in 
[11] show that the collisional energy loss can roughly increase the energy loss by 30-40% 
for RHIC energies. However, it has recently been claimed [12, 13, 14] that for RHIC 
conditions the collisional energy loss may be as important as the radiative one, or even 
dominate at low energies. But an accurate comparison of the two mechanisms of the 
energy losses so far has not been performed, say, in [14] even for the plasma with a 
uniform density the radiative energy loss was calculated incorrectly. The authors have 
used the kinematical suppression factor for the radiative energy loss obtained in [8] which 
strongly overestimates the kinematic suppression [15]. The question of the importance of 
the collisional energy loss is becoming of special current interest in connection with the 
recent data on the non-photonic electrons [16] indicating that the nuclear suppression for 
heavy quarks may be similar to that for light ones. This fact seems to be inconsistent with 
purely radiative energy loss which, at low energies, should be suppressed for heavy quarks 
by the mass effects [17]. Unfortunately, presently for both the radiative and collisional 
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losses the uncertainties in the theoretical predictions are rather big. Say, the results are 
very sensitive to the choice of a s (running or constant [11, 18]), to the infrared effects [11]. 
For clarifying the situation with the relative contributions of the radiative and collisional 
energy loss it is important to perform the calculations with the same parametrizations of 
the coupling constant and the infrared cutoffs. In the present paper we present the results 
of such calculations for expanding QGP for RHIC and LHC conditions. 

To evaluate the radiative energy loss we use the LCPI formalism [4, 6, 7]. It accurately 
treats the mass and finite-size effects, and applies at arbitrary strength of the Landau- 
Pomeranchuk-Migdal (LPM) effect [19]. Other available approaches have limited domains 
of applicability, and can only be used either in the regime of strong (the BDMPS formalism 
[3]) or weak (the GLV formalism [8]) LPM suppression (the GLV approach [8], in addition, 
is restricted to the emission of soft gluons). 

Following Bjorken [2] we evaluate the collisional energy loss for elastic binary collisions. 
However, contrary to the Bjorken analysis we treat accurately the kinematics of the binary 
collisions and fluctuations of the momentum transfer due to the thermal parton motion 
in the QGP. For the infrared momentum cutoff we use the Debye mass included into the 
gluon propagator. In pQCD more accurate treatment of the infrared region of collective 
excitations is possible in the Hard Thermal Loop (HTL) resummation technique [20]. 
However, it is unlikely that the pQCD HTL formalism assuming that g <C 1 is reliable 
for RHIC and LHC conditions when the plasma temperature T < (2 — 3)T C (T c « 170 
MeV is the temperature of the deconfinement phase transition) and g ~ 1.5 — 2. Say, 
in the leading order in g the HTL approximation predicts zero magnetic screening mass. 
However, the lattice calculations [21] show that the magnetic mass may be of the order 
of the electric Debye mass. Also, the HTL approach to the collisional energy loss [22, 23] 
does not incorporate the running a s which turns out to be very important [18]. 

In our calculations of the radiative and collisional energy loss the infrared cutoffs for 
the t-channel gluon exchanges are performed in the same way. This should minimize 
the theoretical uncertainties associated with the collective excitations in the ratio of the 
radiative to collisional energy loss which is of interest in the present paper. Numerical 
calculations are performed with the Debye mass obtained in the lattice calculations. Our 
results show that for RHIC and LHC conditions the effect of the collisional energy loss is 
relatively small, and cannot be crucial for the difference in the nuclear suppression of the 
heavy and light quark jets. For c-quark we find the radiative energy loss which is very 
close to that for the light quarks. For 5-quark the radiative energy loss is significantly 
suppressed at E < 20 GeV, but at high energies it exceeds the charm (and light quark) 
energy loss. Our results show that the observed at RHIC suppression of the non-photonic 
electrons [16] may be naturally described if in the experimentally studied region pt < 8 
GeV the non-photonic electrons are dominated by charm production. Although the pQCD 
calculations [24] predict that the bottom contribution dominates at pr > 5 GeV, this 
possibility cannot be excluded since the pQCD results in the charm mass region are 
very fragile, say, the calculations of Ref. [24] underestimates the normalization of the 
experimental electron spectrum by a factor of about 5 [25]. The same situations occurs 
for the .D-meson production [26]. Thus, presently, it is not clear at all whether the pQCD 
calculations predict correctly the ratio charm/bottom for the kinematical region studied 
at RHIC. 
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2. We begin with the radiative energy loss. For definiteness we consider the case of an 
energetic quark. We assume that the fast quark, Q, produced in a hard process at z — 0, 
passes through a length L of an expanding QGP (we choose the z-axis along the initial 
quark momentum). We define the energy loss as 



AE = E J dxx 



dP 

dx ' 



where E is the initial quark energy, x is the Feynman variable for the radiated gluon, 
and dP/dx is the probability distribution of the induced gluon emission. Since for hard 
gluons with x > 0.5 the jet really does not disappear, from the point of view of the jet 
quenching, a reasonable choice for the upper limit of ^-integration is x max = 0.5. For x m i n 
we use the value m g /E, where hereafter m g is the gluon quasiparticle mass (the quark 
quasiparticle mass will be denoted m q ). 

In the LCPI approach [4] dP/dx can be expressed in terms of the Green's function 
for a two-dimensional Schrodinger equation in the impact parameter space in which the 
longitudinal coordinate z plays the role of time. This Schrodinger equation describes 
evolution of the light-cone wave function of a spurious three-body qqg color singlet system. 
The corresponding Hamiltonian for the medium number density n(z) has the form 

H = ~2M{x)\d P ) ~ % 2 + V (2) 

Here M{x) = Ex(l — x) is the reduced "Schrodinger mass", Lf — 2M(x)/e 2 with e 2 = 
m 2 x 2 + m 2 (l — x) 2 , a 3 (p,x) is the cross section of interaction of the qqg system with a 
medium constituent which reads a 3 (p,x) = \[<J q q(p) + c q q((l — x)p)\ — \o qq -(xp) [27, 4], 
where 

[1 - exp(iq ± p )] 

is the dipole cross section for the color singlet qq pair (Cf,t are the color Casimir for the 
quark and thermal parton (quark or gluon), po is the Debye mass). The gluon spectrum 
may be written as [17, 11] 



v q q(P) = a s C T C F J dq ± (3) 







L 

dx J dx 



dP f , da* H f {x,z) 

dz n{z) — , (4) 



da BH (x,z) f f f 

eff dx ' — = Re J dzi J dz 2 J dpgK v (z 2 , p 2 \z, p)a 3 (p,x)K(z, p\z 1: p x ) =p ^ =Q - ( 5 ) 

z 12 

where K is the Green's function for the Hamiltonian (2), and K v is the Green's function 
for the same Hamiltonian with n = 0, g = a s CV[l + (1 — x) 2 } [2xM 2 (x)]~ 1 • is the 

vertex operator. The dafjj/dx can be viewed as an effective Bethe-Heitler cross section, 
which accounts for the LPM and finite-size effects. It can be represented as [11] 



■togfos) = _ .c Fl i ± (i - ^i Im K l (ftert ) , (6) 

dx ttxM x) J s dp \ ,fp 



ttxM(x) J dp\ Jp 



p=0 
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where the function F is the solution to the radial Schrodinger equation for the azimuthal 
quantum number m — 1 (we omit the argument x) 




Ov 3 (p) 4m 2 -1 J_ 
8Mp 2 L/ 



^,P). (7) 



The boundary condition for F(£,p) reads F(£ — 0,p) = ^/pa3(p)eKi(ep), where K\ is 
the Bessel function. The time variable £ in (7), in terms of the variables z and z\ of 
equation (5), is given by £ = z — Z\\ i.e., contrary to the Schrodinger equation for the 
Green's functions entering (5), (6) represents the spectrum through the solution to the 
Schrodinger equation, which describes evolution of the qqg system back in time. It allows 
one to have a smooth boundary condition, which is convenient for numerical calculations. 

The above equations are written for fixed \xd and a s . For z-dependent po the dipole 
cross section, besides p, depends on z as well. The incorporation of the z-dependence 
of the dipole cross section does not lead to any problems. To account for the effect of 
running a s on the dipole cross section we include a s (q 2 ) in the integrand in (3). The 
inclusion of the running a s in the vertex operator g which corresponds to the emitted 
gluon is a more delicate question since (5) is written in the coordinate representation, 
and there is no explicit dependence on the parton transverse momenta in the vertex 
q — > gq. To generalize the formula (5) to the running a s we use the Schrodinger diffusion 



relation p ~ y (z — zi)/M(x) connecting the longitudinal scale z — z\ and p scale in Eq. 
(5). In terms of the transverse separation between quark and gluon p the transverse 
gluon momentum can be estimated via the uncertainty relation q ~ 1/p. Thus it seems 
quite reasonable to use for the virtuality scale in the splitting vertex the parametrization 
q 2 PS aM(x)/{z — Zi). We adjusted the coefficient a to reproduce in our formalism the 
N = 1 rescattering contribution (which dominates the gluon spectrum) evaluated in the 
diagrammatic approach in momentum representation [15]. For our parametrization of the 
running a s (discussed in Sec. 4) it gives a ps 1.85. 

The induced gluon spectrum can also be used to estimate the effect of the energy 
gain due to absorption of the thermal gluons. For the plasma with fixed (^-independent) 
temperature in the collinear gluon approximation one can obtain for the energy gain 

/dP 
dxx(l + x)n B {Ex) — {y = x/(l + x), E') , (8) 

where E' = E(l + x) is the parton energy after thermal gluon absorption, Ub{p) = 
[e v l T — 1) _1 is the thermal gluon distribution. An accurate evaluation of this effect 
for an expanding plasma in the situation when the plasma density changes strongly at 
the gluon formation scale Lf is a complicated problem. To estimate the effect of gluon 
absorption we have simply used for the thermal distribution the distribution averaged 
over z. Our numerical calculations show that for RHIC and LHC conditions the effect 
of the gluon absorption is suppressed by about two orders of magnitude as compared to 
the gluon emission. For this reason it can be safely neglected, and the difficulties with its 
evaluation are not important for the jet quenching phenomenology. 
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3. Let us now discuss the collisional energy loss. As we said we calculate it through the 
energy transfer in the binary elastic collisions of the fast parton with the thermal quarks 
and gluons treated as free particles. In this approximation the collisional energy loss per 
unit length can be written as 



dE 1 ^ r dp' r dkn p (k) 
~dz~ = 2Ek) J^„ 9p J 2£'(2tt) 3 J 2A;(2vr) 3 



p=q,g 



x 



dk'[l + e p n p (k')} 
2k'(2n) 3 



(2vr) 4 5 4 (P + K-P'- K')u(\M(s,t)\ 2 )9{cu r , 



(9) 



where v xs 1 is the velocity of the fast quark, P = (E, p) and K = (k, k) are the 
momenta for incoming partons, P' = (£", p') and K' = (k', k') are the momenta for 
outgoing partons, u = E — E' is the energy transfer, M(s, t) is the matrix element for 
Qp ~ * Qp scattering (s — (P + K) 2 , t = (P — P') 2 are the Mandelstam variables), 
n q (k) = {e k l T + 1) _1 and n g (k) = {e k l T — 1) _1 are the Fermi and Bose distributions, 
e q = -1, e g = 1, g q = AN c N f , g g = 2(N 2 - 1). In (9) u max is the upper limit of the energy 
loss. Similarly to the radiative energy loss we take u max = E/2. After integrating over 
the p', azimuthal angle of the transverse momentum k^ and k' z , (9) takes the form 



dE 



x 



J 



dz 1QEE'v(2tt) 
dq ± [l + e p n p (k')} oo 
k' J 



p=q,9 



dk z dk±k± n p (k) 
k 



/ 



(\M(s,t)\ 2 )6(u max -u), 



dq z 



K 

t4 + 



p, 



k ' ^rn 2 Q + (P z - q z y + qi 



(10) 



where q = k' — k. The longitudinal component q z is determined from the energy conser- 
vation E + k = E' + k! . At small \t\ the amplitude is dominated by the t-channel gluon 
exchange which gives for the average squared matrix element 



(\M(s,t)\ 2 )^C p 



2™g(|t|) 
(\t\+»l) 2 



(12) 



with C q = N £ N 2 , Cg — 1. The u may be written as (if we take v — 1) 



-t - tk z /E + 2k ± q ± 
2(k - k z ) 



(13) 



Note that in the Bjorken analysis [2] the last two terms in the numerator of (13) have been 
neglected. In this case neglecting the statistical Pauli-blocking and Bose enhancement 
factors one can represent (9) as 



dE 1 /■ 
dz ~ 2(2vr)3 ^„ 9p J 



P=Q,9 



\t\max 

„ nJk) f , . .da 

ik — / dm Tt' 

o 



k 



(14) 
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where \t\ max ~ 2{k — k z )u max . Eq. (14) is convenient for numerical calculations. However, 
at low energies E < 10 GeV, it is not accurate enough. For this reason we use the form 
(10). In numerical calculations we have used accurate formulas for the matrix elements 
[28, 29]. Note that for heavy quark at E <m 2 q /T (in this energy region the heavy quark 
becomes nonrelativistic in the centre of mass system of the binary collision) the value of 
u) max may be smaller than E/2 due to the kinematical limits. It suppresses the energy 
loss at low energies. 

4. To apply our formulas we need to specify the parametrization of a s and mass parame- 
ters. We parametrize a s (Q 2 ) by the one-loop expression which is frozen at some value a{ r 
at Q < Qf r - Previously such a form with a{ r ~ 0.7 was used in the analyses of the low- a; 
structure functions within the dipole approach [30, 31, 27]. A similar parametrization has 
been used in [32] in the analysis of the heavy quark jets. From the analysis of the heavy 
quark energy loss in vacuum the authors of Ref. [32] obtained 



For our parametrization from (15) one can obtain a s (Q < Qf r ) = 0-{ r ~ 0.7, and Qf r ~ 
0.82 GeV (for Aq C d = 0.3 GeV) which agree surprisingly well with the parameters of 
Refs. [30, 31]. In the vacuum the stopping of the growth of a s at Q < Qf r ~ 1 GeV 
may be caused by the nonperturbative effects [32, 33]. In the QGP thermal partons 
can give additional suppression of a s at low momenta [34]. Unfortunately, at present 
there is no robust information on a s (Q 2 ) in the QGP for gluons interacting with the 
energetic (E ^> T) partons which is necessary in our case. Available pQCD calculations 
are performed in the static limit (see, for example, [35, 36, 37] and references therein). 
The running coupling constant obtained in [36, 37] has a pole at Q/Aqcd ~ 3 at T ~ 250 
MeV. Thus, in pQCD, even for the static case, the situation with Q-dependence of the 
in-medium a s is unclear. In the absence of robust analytical theoretical predictions for 
the in-medium a s for fast partons it seems reasonable to estimate a{ r from the lattice 
results on the thermal a s (T). The lattice simulations [38] give a s (T) smoothly decreasing 
from ~ 0.5 at T w 175 MeV to ~ 0.35 at T « 400 MeV. This behaviour of a s in the QGP 
is also consistent with the analysis of the lattice data within the quasiparticle model [39]. 
One can expect that the thermal a s (T) should be somewhat smaller than a{ r . For this 
reason it seems reasonable to use a{ r ~ 0.5 for RHIC and LHC conditions. 

The collisional energy loss of light quarks and gluons is not sensitive to the quark and 
gluon quasiparticle masses. However, this is not the case for the induced gluon radiation 
which is especially sensitive to m g . In the pQCD HTL resummation [20] m q = gT/^/3 

and m g = gT^J (1 + Nf/2)/2. Since the HTL pQCD formulas may be unreliable for RHIC 
and LHC conditions it seems better to to use the results of the lattice simulations. We use 
the quasiparticle masses obtained in Ref. [39] from the analysis of the lattice data within 
the quasiparticle model. For the relevant range of the plasma temperature T ~ (1 — 3)T C 
the analysis [39] gives m q m 0.3 and m g rs 0.4 GeV. 

Besides the quasiparticle masses, we need to specify the the Debye mass which en- 
ters the dipole cross section and the amplitudes of the binary collisions. We perform 
calculations for a fixed and T-dependent jio- In the first case we use the Debye mass ob- 




(15) 
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tained in [39] through the perturbative relation fxo — V2rn g with m g extracted from the 
quasiparticle fit of the lattice data, which gives approximately the T-independent value 
Hr> ~ 0.57 GeV. For the T-dependent parametrization we take the Debye mass obtained 
in the lattice calculations for Nf = 2 [40] which give the ratio fio/T slowly decreasing 
with T (n D /T « 3 at T ~ 1.5T C , fx D /T « 2.4 at T ~ 4T C ). 

For evolution of the QGP produced in AA-collisions we use the Bjorken model [41] 
with the longitudinal expansion which gives the proper time dependence of the plasma 
temperature T 3 r = Tq t (T is the initial plasma temperature). For fast partons produced 
in the central rapidity region of AA-collisions our coordinate z equals the proper time 
r. Thus, we have the number density n(z) oc 1/z for z > r . In the mixed phase the 
fraction of the QGP was calculated according to the 1/r dependence of the entropy density. 
Rescatterings in the hadron phase giving a small contribution have been neglected. As 
in our earlier analysis of the nuclear modification factor [11] we assume that the QGP is 
in the thermal and chemical equilibrium (we take Nf = 2.5). For RHIC we performed 
calculations for T = 297 MeV and r = 0.5 fm which correspond to the initial conditions 
in successful hydrodynamic description of Au + Au collisions at y/s = 200 GeV [42], and 
agree with the total entropy extracted from the charged particle multiplicity. For Pb + Pb 
collisions at LHC for ^/s = 5500 GeV we use T = 350 MeV (with the same r = 0.5 
fm), which was obtained from the extrapolation of the RHIC data on the charged particle 
distribution to the LCH energy performed in [43]. 

5. We present the numerical results for a{ r = 0.5. In Figs. 1, 2 we plot the radiative and 
collisional energy losses for the light quark and gluon for RHIC and LHC conditions. To 
illustrate the effect of the running coupling constant we present the results for the case 
a s(Q) — ttf r as well. The higher panels correspond to the T-independent he, = 0.57 GeV 
[39], and the lower panels to the T-dependent Debye mass from Ref. [40]. The calculations 
were performed for L = 5 fm which is the typical parton pathlength in the QGP (and 
mixed) phase with life-time about Ra ~ 6 fm for the central heavy ion collisions. One sees 
that the fraction of the collisional energy loss is relatively small. At E ~ 10 GeV the ratio 
AE co i/AE rac i is about 0.3-0.4 for quarks and 0.2-0.3 for gluons. The smaller fraction of 
the collisional energy loss for gluons results from the additional color factor Ca/Cf = 9/4 
for the radiative energy loss for gluons (which is absent for the collisional energy loss). For 
this reason the effect of the collisional energy loss on the nuclear modification factor should 
be weaker in the kinematical regions where the high-px hadron spectra are dominated by 
gluon jets. The fraction of the collisional energy loss drops as energy increases. The 
T-dependent parametrization of the Debye mass gives somewhat smaller AE. It is due 
to suppression of the rescatterings in the initial high-temperature region with z <j 1 fm 
where \iv may be about 1 GeV. However, one can see that the sensitivity of the results 
to the Debye mass is relatively weak. The results shown in Figs. 1, 2 demonstrate clearly 
the importance of the running a s for both the radiative and collisional energy loss. It 
leads to flattening the energy losses at high energies. 

In Fig. 3 we show the results for the charm (m c =1.2 GeV) and bottom (m& = 4.5 
GeV) quarks obtained with running a s for the Debye mass from [40]. The difference in 
the radiative energy loss for light and charm quarks is small. The charm radiative energy 
loss is only suppressed by ~ 10% at E ~ 10 GeV compared to the light quarks. For 
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the bottom quark the mass suppression at low energies is significant. Note that at high 
energies the radiative energy loss for the bottom quark becomes larger than that for the 
charm quark, and the charm contribution exceeds slightly the light quark one. This fact 
is connected with emission of gluons at moderate values of x where the induced radiation 
for heavy quarks turns out to be enhanced in the so-called diffusion regime when L <C Lf 
[44] . A detailed discussion of this effect reflecting a complicated interplay of the finite-size 
and mass effects in the induced gluon emission will be given elsewhere. 

A small difference in the energy loss for the light and charm quarks show that one 
can expect approximately the same nuclear suppression for light and c-quarks. This is 
in contradiction with the considerable suppression of the induced gluon radiation from 
charm predicted in [45]. However, in [45] there was not performed any accurate evalua- 
tion of the mass effects in the induced radiation. To obtain the heavy quark spectrum 
the authors multiplied the BDMPS spectrum obtained for massless partons in the oscil- 
lator approximation by the suppression factor defined as the ratio of the vacuum gluon 
emission spectra for heavy and light quarks which evidently has noting to do with the 
mass modification of the induced gluon radiation. The enhancement of the energy loss 
for heavy quarks at high energies is absent in [45]. 

To study the infrared sensitivity of our results we also performed computations for 
m g = 0.2 and m g = 0.6 GeV. These values of m g give reasonable the lower and upper 
limits of the infrared cutoff for the induced gluon emission for RHIC and LHC conditions 1 . 
At E ~ 10 GeV for m g = 0.2 GeV AE rad is bigger by ~ 20 - 30% than that for m g = 0.4 
GeV, for m g = 0.6 GeV the effect is of opposite sign. The effect of variation of the gluon 
mass in the above interval of m g becomes small at higher energies (< 10 — 15% at E > 40 
GeV). 

To study the effect of variation of a s we have also performed calculations using for 
a{ r the values 0.7 and 0.35. The first one neglects the in-medium suppression of coupling 
constant, and the second one, in the light of the lattice results, can be viewed as an low 
bound for a{ r for RHIC and LHC conditions. Numerically we obtained approximately 
the same fraction of the collisional energy loss as for a{ r = 0.5. Thus, for the reasonable 
bandwidths in m g and a{ r the fraction of the collisional energy loss is small. 

Note that due to the dominance of the radiative energy loss modeling the jet quenching 
with the collisional energy loss alone within the model of a particle undergoing Brownian 
motion described by the Fokker-Planck equation [12] does not make sense. Evidently 
for an accurate evaluation of the nuclear modification factor the radiative and collisional 
energy losses must be treated on an even footing. One can expect a nontrivial interplay of 
these two effects, say, it is clear that they cannot be additive, since the collisional energy 
loss will suppress the effective in-medium gluon formation length. The quantum nonlocal 
character of the induced gluon radiation makes the problem very complicated even at the 

1 Note that the analysis of the low- a; proton structure function within the dipole BFKL equation 
[27, 31] gives the value of the effective gluon mass for gluon emission in the parton-nucleon interaction 
about 0.75 GeV. This value agrees well with the natural infrared cutoff for gluon emission in the vacuum 
m g ~ 1/R C , where R c w 0.27 fm is the gluon correlation radius in the QCD vacuum [46]. One can expect 
that in AA-collisions the infrared cutoff will be approximately the same only for gluon emission in the 
developed mixed phase and for fast gluons with Lf >• L which give relatively small contribution to the 
total radiative energy loss. 
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level of the radiative energy loss. Presently the distribution in the induced energy loss 
which is necessary for evaluation of the nuclear modification factor is usually calculated 
assuming the independent gluon emission [47] which, however, has not any serious the- 
oretical justification. Since the fraction of the collisional energy loss is small, it seems 
reasonable to treat it as a perturbation. At qualitative level, neglecting the nonadditivity, 
one can incorporate the collisional energy loss into this model by a small renormalization 
of the QGP density according to the change in the AE due to the collisional energy loss. 
In [11] we have described reasonably well the RHIC data on the nuclear modification 
factor in Au + Au collisions at y/s = 200 GeV by the induced gluon radiation alone with 
a{ r = 0.7. Inclusion of the collisional energy loss will require somewhat smaller value of 
a{ r . The results of this analysis will be presented elsewhere. 

Note that although our approach does not treat accurately the region of soft momen- 
tum transfer q < we can expect that this inaccuracy should be small. Indeed, say, 
the relative contributions into collisional energy loss of the soft region with q < 2/ip eval- 
uated in our approach and in the HTL pQCD approach [23] with accurate treatment of 
the collective excitations are close. Also, even for low parton energy E ~ 5 — 10 GeV 
the soft region gives relatively small effect (about 30% for T ~ 250 MeV). In any case, 
since the infrared effects should modify the dipole cross section (3), which controls the 
induced radiation, and the probability of the collisional energy transfer approximately in 
the same way one can expect a good stability of the ratio of the two mechanisms against 
the inaccuracy in the soft momentum region. 

6. In summary, we have performed the comparison of the radiative and collisional energy 
losses of energetic quarks and gluons in an expanding quark-gluon plasma for RHIC and 
LHC energies. The radiative energy loss has been calculated within the LCPI approach 
[4]. To evaluate the collisional energy loss we have used the Bjorken model of elastic 
binary collisions with an accurate treatment of kinematics of the binary collisions. The 
calculations have been performed with the same infrared cutoffs and parametrization of 
the coupling constant for the radiative and collisional energy loss, which is important 
for minimizing the theoretical uncertainties in the ratio of the radiative and collisional 
contributions. Our numerical results demonstrate that for RHIC and LHC conditions the 
fraction of the collisional energy loss is relatively small, and decreases with energy. For 
gluons it is smaller than for quarks. 

Our calculations show that the difference in the radiative energy loss for charm and 
light quarks is small. For this reason the nuclear modification factor for light hadrons and 
D-mesons should be approximately the same in the kinematical region where the light 
hadron spectra are dominated by the quark jets. At sufficiently large energies the heavy 
quark energy loss becomes bigger than that for light quark. 
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Figures 
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Figure 1: The light quark radiative (solid line) and collisional (dashed line) energy loss 
for RHIC (left) at ^ = 200 GeV and LHC (right) at y/s = 5.5 TeV conditions for L = 5 
fm. The thick curves correspond to the running a s , and thin curves to a s = 0.5. The 
higher panels show the results for the T-independent Debye mass nr> ~ 0.57 GeV [39], 
and the lower panels for the T-dependent Debye mass from the lattice calculations [40]. 
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Figure 2: The same as in Fig. 1 but for gluon. 
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Figure 3: The charm (thick curves) and bottom (thin curves) quark radiative (solid line) 
and collisional (dashed line) energy loss for RHIC (left) at y/s = 200 GeV and LHC (right) 
at y/s = 5.5 TeV conditions for L = 5 fm, m c = 1.2 GeV, m b = 4.5 GeV. The dotted line 
shows the radiative energy loss for light quark. The calculations were performed with the 
running a s and the T-dependent Debye mass from [40]. 
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